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A method of characterizing fluorescent molecules or other particles using generating functions 

This invention relates to the field of fluorescence spectroscopy, and more particularly to a 
method for determining characteristic physical quantities of fluorescent molecules or other 
particles present in a sample. 

The primary data of an experiment in fluorescence correlation spectroscopy (FCS) is a 
sequence of photon counts detected from a microscopic measurement volume. An essential 
attribute of the fluorescence correlation analysis is the calculation of the second order 
autocorrelation function of photon detection. This is a way how a stochastic function (of 
photon counts) is transformed into a statistical function having an expected shape, serving as a 
means to estimate some parameters of the sample. However, the calculation of the 
autocorrelation function is not the only way for extracting information about the sample from 
the sequence of photon counts: Further approaches are based on moment analysis and analysis 
of the distribution of the number of photon counts per given time interval (Qian and Elson, 
Pioc. Natl. Acad. ScL USA, 87: 5479 - 483, 1990; Qian and Elson, Biophys. J. 57: 375 - 380, 
1990). 



The intensity of fluorescence detected from a particle within a sample is not uniform but 
depends on the coordinates of the particle with respect to the focus of the optical system. 
Therefore, a reliable interpretation of measurements should account for the geometry of the 
illuminated measurement volume. Even though the calculation of a theoretical distribution of 
the number of photon counts is more complex for a bell-shaped profile than for a rectangular 
one, the distribution of the number of photon counts sensitively depends on values of the 
concentration and the specific brightness of fluorescent species, and therefore, the measured 
distributions of the number of photon counts can be used for sample analysis. The term 
„specific brightness'* denotes the mean count rate of the detector fiom light emitted by a 
particle of given species situated in a certain point in the sample, conventionally in the point 
where the value of the spatial brightness profile function is unity. 



The first realization of this kind of analysis was demonstrated onthe basis of moments of the 
photon count number distribution (Qian and Elson, Proc. Natl. Acad. Sci. USA, 87: 5479 - 



483, 1 990); The fc-th factorial moment of the photon count number distribution P(rt) is defined 

as ■ v i ■- ■ . ■ : , . • 

» 1 *.''*, - - ^ J ' ' 

: * * J ."'.*, l , ' * 

In turn, factorial moments are closely related to factorial cumulahts, •• *. < " 

Or ■ - ■ • • * ' - . " : ., ■ ; ' _ ' ■ '* • 

X^K-ItCfr'K^F,. . ; ''V, . ; (3) ■•■ 

( Cf s are binomial coefficients, and AT^s are cumulants.) The basic expression used in moment 
analysi$> derived for ideal solutions, relate A-th order cumulant to concentrations (c/) and 
specific brightness values ) ' 

' ' • . • ' . ■ - ■ ■ 

Here, Xk l% monient of the relative spatial brightness profile B{r): ■ 

x, = jB k (r)dv. r ' j - ' (5) 

, • • . ■ • ' ■ ' ' ■ • . - - - ■■ . 

** . ■'»■■* j 

* ^* . * 

* * / " ' 

Usually in FCS, the unit of volume and the unit of B are selected -which yield Xi - Xi'—'l- 
After selecting this convention, concentrations in the equations are dimensionless, expressing 
the mean number of particles per measurement volume, and the specific brightness of any 
species equals the mean count rate from a particle if situated in the focus divided by the 
numeric value of 5(0). The value of this constant is a characteristic of optical equipment. It 
can be calculated from estimated parameters of the spatial intensity profile (see below). Qian 
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and Elson used experimental values of the first three cumulants to determine unknown 
parameters of the sample. The number of caimulants which can be reliably determined from 
experiments is usually three to four. This sets a limit to the applicability of the moment 
analysis. 

The idea behind the so-called fluorescence intensity distribution analysis (FIDA), which has in 
detail been described in the international patent application PCT/EP 97/05619 (international 
publication number WO 98/16814), can be well understood by imagining an ideal case when a 
measurement volume is uniformly illuminated and when there is almost never more than a 
single particle illuminated at a time, similar to the ideal situation in cell sorters. Under these 
circumstances, each time when a particle enters the measurement volume, fluorescence 
intensity jumps to a value corresponding to the brightness of a given type of particles. 
Naturally, the probability that this intensity occurs at an arbitrary time moment equals the 
product of the concentration of a given species and the size of the measurement volume. 
Another fluorescent species which may be present in the sample solution produces intensity 
jumps to another value characteristic of this other species. In summary, the distribution of 
light intensity is in a straightforward way determined by the values of concentration and 
specific brightness of each fluorescent species in the sample solution. 

It is assumed that the light intensity reaching the detector from a particle as a function of 
coordinates of the particle is constant over the whole measurement volume, and zero outside 
it. Also, it is assumed that the diffusion of a fluorescent particle is negligible during die 
counting interval 7*. In this case, the distribution of the number of photon counts emitted by a 
single fluorescent species can be analytically expressed as double Poissonian: the distribution 
of the number of panicles of given species within this volume is Poissoniaai, and the 
conditional probability of the number of detected photons corresponding to a given number of. 
particles is also Poissoniah. The double Poissonian distribution has two parameters: the mean 
number of particles in the measurement volume, c and the mean number of photons emitted 
by a single particle per dwell time, qT. The distribution of the number of photon counts n 
corresponding to a single species is expressed as 

Ptoc.q)*±£<«&ZL'r-*-. ; ( 6) 



where m runs over the number of molecules in the measurement volume. Jf-P\(n) denotes the 
distribution of the number of photon counts from species /, then the resultant distribution P(n) 
is expressed as - / ; • •.. v. 



This me^ that /'(w) can be calculated as a convolution of tiie sdies of distributions Pi{n). 



Like in FCS, the rectangular sample profile is a theoretical model which can hardly be applied 
in experiments. One may divide the measurement volume into a great number of volume 
elements and assume that within each of them, the intensity of a molecule is constant 
Contribution to photon count number distribution fn>m a volume element is therefore double 
Poissonian with parameters cdV snd qfB(r). (Here q denotes count rate from a molecule in a 
selected standard position where B~ 1, and B(r) is the spatial brightness profile function of 
coordinates.) The overall distribution of the number of photon counts can be expressed as a 
convolution integral over double Poissonian distributions. Integration is a one-dimensional 
rather than a three-dimensional problem here, because the result of integration does not 
depend on actual positions of volume elements in respect to each other. Figuratively, one may 
rearrange the thite-dimensional array of volume elements into a one-dimensional array, for 
example in the decreasing order of the value of ifc " 



. In a number of first experiments described . in the international patent application PCT/EP 
97/05619, the photon count number distribution was indeed fitted, using the convolution 
technique. The sample model consisted of twenty spatial sections, each characterized by its 
volume V] and brightness By However, the t^hnique described in this patent application is 

7 slow and inconvenient in cases involving a high number of samples to be analyzed, like in 
diagnostics or drug discovery, or in analyzing distribution functions involving more than a 
single argument ■„.,' , •• ' r " 



Therefore, it is ah object of the present invention to present a convenient and much fester 
technique for analyzing fluorescence intensity fluctuations. 




According to the present invention there is provided a method for characterizing fluorescent 
molecules or other particles in samples, the method comprising the steps of: 

a) monitoring fluctuating intensity of fluorescence emitted by the molecules or other particles 
in at lea$t one measurement volume of a non-uniform spatial brightness profile by 
measuring numbers of photon counts in primary time intervals by a single or more photon 
detectors, 

b) determining at least one distribution function of numbers of photon counts, ^(n), from 
the measured numbers of photon counts, 

c) determining physical quantities characteristic to said particles by fitting the experimentally 
determined distribution function of numbers of photon counts, wherein the fitting 
procedure involves calculation of a theoretical distribution function of the number of 
photon counts P(n) through its generating function, defined as = ^(n). 



n 



The formal definition of the generating function of a distribution f(n) is as follows: 



O{0 = £4" W (8) 

if=0 w 



What makes the generating function attractive in count number distribution analysis is the 
additivity of its logarithm: logarithms of generating functions of photon count number 
distributions of independent sources, like different volume elements as well as different 
species, are simply added for the calculation of the generating function of the combined 
distribution because the transformation (8) maps distribution convolutions into the products of 
the corresponding generating functions. 

In a particular preferred embodiment, one might monitor the fluctuating fluorescence intensity 
in consecutive primary time intervals of equal width Typical primary time intervals have a 
width in the order of several tenth of microseconds. The total data collection time is usually 
several tenth of seconds. 

In a further preferred embodiment, numbers of photon counts {n-,} subject to determination of 
a distribution function in step b) are derived from numbers of photon counts in primary 



time intervals (Nj) by addition of numbers of photon counts from primary time intervals 
according to a predetermined rule. One might e.g. be interested in choosing numbers of 
photon counts {ni} subject to detemnnation of a distribution function P(*) which are 
calculated from the numbers of photon counts in primary time intervals {N,} according to the 

rule n; ^TN^, where M is an integer number expressing how many times , the time 
interval in which {m} is deteimined is longer than the pimary time interval. 

In a further embodiment, numbers of photon counts {n*} are derived from predetermined 
; primary time intervals according to a rule in which primary time intervals are separated by a 

time delay. In particular^ the following rule can be applied: n, = ^ (N^^ + N ), where 

M and L are positive integer numbers, {n\} are numbers of photon counts subject to 
determination of a distribution function Rn) , and {Nj} are the numbers of photon counts in 
• primary time intervals. 

In some cases, it might be preferred to determine not only a single distribution function in step 
b), but rather a set. of distribution functions P(nj. These can be detennined according to a set 
of different rules, said set of distribution functions being fitted jointly in step c). As an 
example a set of distribution functions with different values of M and/or L might be fitted 
'■ jointly.' ..' - '* 

Typical physical quantities which might be detennined in step c) according to the present 
invention are concentration, specific brightness and/or diffusion coefficient of molecules or 
other particles. ; 

■ ■ v ~ . l i * * . ■ . i - - ■ f 

; - ' * "» * * , " " 

In a further preferred embodiment, the generating function is calculated using the expression 
Gfe) » exp[ \dqc{q) jd 3 r(e ii ~ i) '' m * - 1)], where c(g) is the density of particles with specific . 

> 1 ' - 

brightness q, T is the length of the counting interval, and B(r) is the spatial brightness profile 
as a function of coordinates. 




Applying the definition (8) to formula (6) with c-» cdV and q qB(r) , the contribution 
from a particular species and a selected volume element </F can be written as 

, T ' t . ' - r f+ * 

' 1 ' " * ' ^ -r ' 1 ^ " " 

Therefore, the generating function of the total photon count number distribution can be 
expressed in a closed form 

G(l) = exp[EcJ( e ( ^ raw -lW- ■ (10) 

In * ■ ' ■ " ' ■ 

Numeric integration according to Eq. (9) followed by a fast Fourier transform is the most 
effective means of calculating - the theoretical -distribution P(n) corresponding to a given 
sample (i.e., given concentrations and specific brightness values of fluorescent species). If one 
selects then the distribution P(n) and its generating function G(q>) are interrelated by 

the Fourier transform. Therefore, it is particularly preferred to select the argument of the 
generating function in the form $= e V and to use a fast Fourier transform algorithm in 

calculation of the theoretical distribution of the number of photon counts out of its generating 
function. 

When calculating the theoretical . distribution P(n) in step c) according to the present 
invention, the spatial brightness profile might be modeled by a mathematical relationship 
between volume and spatial brightness. In particular, one might model the spatial brightness 

profile by the following expression: ~ = a x x + c 2 * J + a 3 x 3 , where dV denotes a volume 

element, x denotes logarithm of the relative spatial brightness, and ai, a 2 and ai are empirically 
estimated parameters. 

Some fluorescent species may have a significantly wide distribution of specific brightness For 
example vesicles, which are likely to have a significantly broad size distribution and a random 
number of receptors, may have trapped a random number of labeled ligand molecules. In order 
to fit count number distributions for samples containing such kind of species, it is useful to 





MUENCHEN 05 



23-12-99 :■ 



llilllilllll 



• 



34-297- 



+49 89,2 



modify Eq. (10) in the following manner. The assumption is made that the distribution of 
brightness of particles q within a species is mathematically expressed as follows: 



(11) 



This expression has been selected for the sake of convenience: all moments of this distribution 
be analytically calculated, using the following formula: 



can 



CO 



x e 



Tja + 1) 



(12) 



It is straightforward to derive the modified generating function of a photon count number 
distribution. One can rewrite Eq. (9) as follows: . 1 "■ .- 




TO 

G(|) = exp[£ c, jdV jdqp(q,a, ,b, )(« 



where 



(13) 



piq;a,b} = 



TO 



-q e T 



(14) 



The integral over q can be performed analytically: 



i 



-1 



1' 



(15) 



The parameters a, and i, are related to the mean bri 
distribution aj by 



2 » 



6. = "IT 



(16) 



In the range of obtained count numbers, the probability to obtain a particular count number 
usually varies by many orders of magnitude, see for example the distribution of Fig, 1. 
Consequently, the variance of the number of events with a given count number has a strong 
dependence on the count number. To determine weights for least squares fitting, one may 
assume that light intensities in all counting intervals are independent Under this assumption, 
one has a problem of distributing M events over choices of different count numbers n 9 each 
particular outcome having a given probability of realization, P(n). Covariance matrix elements 
of the distribution can be expressed as follows: 

(17> 

where M is the number of counting intervals per experiment. 

For a further simplification, one may ignore the second term on the right side of Eq. (17), 
which can be interpreted as a consequence of normalization. In this case, the weights simply 
equal to the inverse values of the diagonal covariance matrix elements 



Dispersion matrix (17) corresponds to the multinomial distribution of statistical realizations of 
histograms. The Poissonian distribution, with the constraint that the total number of counting 
intervals U is fixed, will lead to the multinomial distribution. This is the rationale behind 
using Poissonian weights as given in Eq, (1 8). 

Let i\ be the expectation value of the number of events of counting * photons and let 
N - ]T«* be their sum. Let m be a statistical realization with M - Y m k . Assume that 
realizations m^n^Jgi obey Poissonian statistics 

* * 

^ 

%^,K) Jyr, «•», (19) 
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where we have introduced the • notation^ 4 n^'iC 

! w, !K =[m 0 , m, ,K ) . The probability of having the total of M events is 



and 



P{M) = 



N! 



(20) 



The conditional probability of having ^,/m,,K events if tiiere is a total of M events is 



Pim^^^M)- 



or 



[p*>P\& f 



(21) 



This is the multinomial distribution where we have introduced: 



Pk = 



M 



and C ^ * M are multinomial coefficients. 



lii general, a linear, or linearized : least squares fitting returns not only the values of the 
estimated parameters* but also their covariance matrix, provided the weights have been 
meaningfully set: It may turn out to be possible to express ..the statistical errors of the estimated 
parameters analytically in some simple cases (e.g., for the rectangular sample profile and 
single species) but in applications at least two-component analysis is usually of interest. 
Therefore, one may be satisfied with the numerical calculations of statistical errors. In 
addition to the „theoretical" errors with the assumption of non-correlated measurements (Eq. 
(17)), in some cases statistical errors have been estimated empirically, making a series of 
about a. hundred FIDA experiments on identical conditions. As a rule, empirical errors are 
higher than theoretical ones by a factor of three to four. Empirical errors appear to be closer to 
the theoretical ones in scanning experiments. Therefore we are convinced that the main reason 



of the underestimation of theoretical errors is the assumption of non-correlated measurements. 
Table 1 compares statistical errors of parameters estimated by fitting a photon count number 
distribution (FID A) according to the present invention and by the moment analysis. Error 
values are determined through processing a series of simulated distributions. The present 
invention is ovenvheimingly better than die moment analysis if the number of estimated 
parameters is higher than three. 



Table 1 







Primary 


Number 


Number of 


' Parameter 


Value 


Percent 


Percent I 


Data 


time 


of species 

; 


e$timated 


specification 


fas in 


error of 


error of 


collection 


interval, 




^ '*' \ 

parameters 




kHz) 


FIDA 


moment 


time, s 


Ms 










(inventi- 


analysis 














on) 




JO.O 


40.0 


l 


2 


c 


0.5 


0.59 


0.54 












60.0 


0.56 


0.51 


10.0 


40.0 


2 


3 




0.05 


2.00 


2.71 










9i 


150.0 


1.54 


1.89 












3.0 


0.53 


0.62 










q2 


5.0 (fixed) 






10.0 


40.0 


2 


4 


c * 


0,05 


2.26 


4.99. . 


. ■ ' 








<li 


150.0 


1.63 


2.53 


i '* , 








c 2 


3.0 


3.18 


17.8 












5.0 

— - — , i. 


3.35 


14.9 



Confocal techniques may be applied to a wide field of applications, such as' biomedicine, 
diagnostics, high through-put drug screening, sorting processes such as sorting of particles 
like beads, vesicles, cells, bacteria, viruses etc.. The conjugate focal (confocal) technique is 
based on using a point source of light sharply focused to a diffraction-limited spot on the 
sample. The emitted light is viewed through a spatial filter (pinhole) that isolates the viewing 
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area to that exactly coincident with the iUuminating spot Thus, the illumination and detection 
apertures are optically conjugated with each other. "Light originating from focal planes other 
than that of the objective lens is rejected, which effectively provides a very small depth of 
field therefore, in a particular preferred embodiment of the present invention, in step a) a 

T r ' 4..'. .' ' ' . ' £ + * 

confbeal microscope is used for mdrritorkg the intensity of fluorescence. In order to achieve a 
high signal-to-noise ratio, it is useful to monitor the intensity/ of fluorescence using an 
apparatus that comprises: a radiation source (12) for providing excitation radiation (14), an 
objective (22) for focussing the excitation radiation (14) into a ^ measurement volume (26), a 
detector (42) for detecting emission radiation (30) that stems 6pm the measurement volume 
(26), and an opaque means (44) positioned in the pathway (32) bf the. emssion radiation (30) 
or excitation radiation (14) for erasing the central' part of the emission radiation (30) or 
excitation radiation (14). It mightbe particularly preferred to use an optical set-up described in 

detail in Figure 9. 

In a further preferred embodiment, the method according to the present invention is applied to 
,fit a joint distribution of photon count numbers. In experiments, fluorescence from a 
microscopic volume with a fluctuating number of molecules is monitored using an optical set- 
up (e g. a confocal. microscope) with two detectors. The two detectors may have different, 
polarizational or spectral response. In one embodiment, concentrations of fluorescent species 
together with two specific brightness values per each species are determined.. The two : 
dimensional fluorescence intensity distribution analysis (2D-FIDA) if used with a polarization 
cube is a tool which can distinguish fluorescent species with different specific polarization 
ratios. This is a typical example of a joint analysis of two physical characteristics of single 
molecules or other particles, granting a significantly improved reliability compared to 
methods focussed on a single physical characteristic. 

*■' ' t m. . ' t ' 



.1 



1 . ■ - - - J ■ f 

In order to express the expected two-dimensional distribution of the number of photon counts, 
it is favourable to use the following assumptions: (A) Coordinates of particles are random and 
independent of each other. (B) Contribution to fluorescence intensity from a particle can. be- 
expressed as a product of a specific brightness of the particle and a spatial brightness profile 
function characteristic to the optical equipment (C) A short counting time interval T is 



selected, during which brightness of fluorescent particles does not significantly change due to 
diffusion. 



At first, a joint distribution of count numbers fiom a single fluorescent species and a single 
small open volume element dV is expressed. The volume element is characterized by 
coordinates r and spatial brightness B(r), and the fluorescent species is characterized by its 
specific brightness values <?, and By 9l and * 2> mean photon count rates by two 
detectors from a particle situated at a point where B(r) = 1 are denoted A convenient choice is 
to select a unit of 5, as usual in FCS, by the equation = x% > where Xi = J**(r)rf 3 r. If the 
volume element happens to contain m particles, then the expected photon count numbers per 
time interval T from the volume element are mq x TB(x) and mq,TB{x) , while the distribution 

of numbers of photon counts from m particles P(n„n 2 \'m) is Poissonian for both detectors 
independently: 



Fiom the other side, under assumption (A), the distribution of the number of particles of gi 
species in the volume element is Poissonian with mean cdV , c denoting concentration: 



P^JA^' (23) 



The overall distribution of the number of photon counts from the volume element can be 
expressed using Eqs.22 and 23: 



> r " •■ 
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(24) 



As in the on^dimensional case described above, a useful representation of a distribution of 

. ■ • '■ . ■ i ■ ■ ' , j. ' 

* . - ' ■ ' ' ' ' " . " ' t • * ' ' r, * f •■ . ' . ' 

numbers of photon counts />(£»*) is its generating function, defined as 



(25) 



The generating 




as 



n 2 ! 



(26) 



In particular, if one selects = exp(i<p, ), then the distribution ^ and its generating 
function £?(?,, <Pr) are interrelated by a 2-dimensional Fourier transforai. What makes the 

" * . . ' ' * n ■ T T ■ i \ 

j* ' ■* * - * "' * "i - 

generating function attractive in photon counrnumber distribution analysis is the additivity of 
its logarithm: logarithms of generating functions of photon count number distributions of 
independent sources, like different volume elements as welt as different species, are simply 

istribution. Therefore, the generating function of 
the overall distribution of the number of photon counts can be expressed in a closed form: 




15 

G&^^expKSr-W*^^^ (27) 

In this formula, a contribution from background count rates, Xi by detector 1 and Xz by 
detector 2, as well as contributions from different fluorescent species, denoted by the subscript 
i have been integrated. Numeric integration according to Eq. 27 followed by a fast Fourier 
transform is a very efficient means for calculation of the theoretical distribution P^n^) 

corresponding to a given sample (i.e. given concentrations and specific brightness values of 
fluorescent species). 

The spatial brightness function is accounted through the spatial integration on the right side of 
Eq. 27. The three-dimensional integration can be reduced to a one-dimensiorial one by 
replacing three-dimensional coordinates r by a one-dimensional variable, a monotonic 
function of the spatial brightness JJ(r). A convenient choice of the variable is 
x = ln[5(0) / B(r)]. A sufficiently flexible model of the spatial brightness profile is presented 
by the following expression: 

dV 

— <**(! + a,x + tf 2 x 2 ). (28) 

In the interval of obtained count numbers, the probability to obtain a particular pair of count 
numbers usually varies by many orders of magnitude. Consequently, the variance of the 
experimental distribution has also a strong dependence on the count numbers. To determine 
weights for least squares fitting, for simplification it is assumed that coordinates of particles in 
all counting intervals are randomly selected. (This means one ignores correlations of the 



-; 1 



coordinates in consecutive costing inters.): Under this assumption, one has a problem of 
distributing M events over choices of different pairs of count numbers m, * each particular 

outcome h^^^ 

distribution can be expressed as follows: „ • • 



. P(n. . ».15 f n, ,n', )5 ,* a ) - ffa ^("'i . ""'» ) (29) 
{AP(n.^)AP(nVi^2))= 7 ; M ... 



where Af is the number of counting intervals per experiment. 



For a further simplification, one may ignore the second term on the right side of Eq.(20) s 
which can be interpreted as a consequence of normalization. In this case, the weights simply 
equal to the inverse values of the diagonal covariance matrix elements 



^ )= y^H)' : - : (30) 



In general, a linearized least squares fitting algorithm returns not only values of estimated 
parameters, but also their covariance matrix, provided weights have been meaningfully set. In 
addition to «'theoretic a l'Verrors:correspondmg to the assumption of uncorrelated measurements 
(Eq. (30)), in some cases statistical errors have been determined empirically, making a series 
of about 100 2D-FIDA experiments identical conditions. As a role, empirical errors are 
higher than theoretical ones by a factor of three to four . The main reason of underestimation of 
theoretical errors is most likely the assumption of uncorrelated measurements. 



Even though the assumption of uncorrelated measurements yields underestimated error 
values, it is a very useful theoretical approximation, allowing to compare accuracy of analysis 
under different experimental conditions as well as different methods of analysis. Also, this 
approximation provides an easy method of data simulation which is a useful tool in general, A 
simple and very fast method of data simulation is calculation of the expected event number as 
a function of photon count numbers and addition of random Poisson noise to each event 
number independently. 

In Table 2 theoretical errors of one-dimensional and two-dimensional fluorescence intensity 
distribution analysis according to the present invention are presented in two selected cases of 
two fluorescent species. In both cases the ratio of specific brightness values of the two species 
is three. In the case of 2D-FIDA, it is assumed that spectral sensitivities of the two detectors 
are tuned to different species. In both cases data collection time of 10 s, time window of 40 ^s 
and background count rate of. 1 kHz are assumed. Note that the statistical errors of the 
estimated parameters are significantly lower in the 2D-FDDA example. 





foENCHENi 05 



.3-12-99 : I 



4297-> ' +4-9 89 23994; 



18 



Table! 



Method Parameter specification 


Parameter value (selected) 


Percent, error 


1D-HDA dV 




6.6 


(according 02V 


0.5 ' / 


4.9 


, . .. ' - ' ' . ■ V ' 

to , the qi 


60 kHz 


' 2.2 ' ' 


* * ■ . ■;• 

invention) qz 


20 kHz 


10,2 


2D-FIDA CiV 


0.5 


1.1 


(according C2V 




1.1 ' 


■■■■■■■ ■ i • • ' • ■'■ •' V--"- 

" V - , * ' ' ' ■ 

to the q*i 


60 kHz 


0.77 

L - 


* ' ■ 1 - , * * ' '* ' 

invention) qA2 


20 kHz 


1.2 


qsi 


20 kHz 


1.2 


qB2 


60 kHz 


0.77 - ; 




The two-dimensional fluorescent intensity distribution analysis according to the present 
invention is in the following compared to the prior art moment analysis. 

Factorial moments of the distribution P{n\/ti) are defined as 



i 



(«, -*)!<«, * 0! 



(31) 



Factorial moments are related to factorial cumulants K k} 



(32) 



IBiiiiSiii^ffll 



C denotes binomial coefficients. Cumulants can be expressed through concentrations and 
specific brightness values by a simple relation 



(33) 




If the unit of B is selected by the equation %\ - Xi » Xi has the meaning of the sample volume, 
denoted by V, and Eq. 33 can be Yfritten as 



(34) 



where y denotes a series of constants characterizing the brightness profile: 




(35) 



Tne principle of moment analysis is to determine values of a few cumulants from an 
experiment and solve a system of Eqs. 33 in respect to unknown concentrations and brightness 
values. 



Li Table 3 statistical errors of 2D-FIDA (present invention) and 2D-MAFED (prior art) are 
presented determined by generating a series of 30 random distributions of count numbers, 
simulated for identical "samples", thereafter applying 2D-FIDA and 2D-MAFID, and 
detennining the variance of estimated parameters in both cases. Note a tendency that the 
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advantages of 2D-FIDA compared to 2D-MAFBD increase with the number of parameters to 
be estimated. 



Tabled 



\ 



Data 
collec- 
tion . 
time, s 



Number 

of 
species 



Number of Specifica- 
estimated tioxi of 
parameters cumulants 

used in 
MAFID 



1 



Specifi- 

cation of 

■ ■ 'i. > 

parame- 
ters 



True 
values 



Percent 
error of 
F3DA 



Koi 
Kio 
Ku 



cV 

qA 

4b 



0.5 

60 kHz 
40 kHz 



1.1 
1.2 
1.3 



Percent 
error of 
MAFID 



1.3 
L5 
1.6 



id 


2 


A 

4 


; . 




0.5 


1,1 : 


2.5 










C2V- 


0 5 


1.0 






. * 






QAl 


60 kHz 


0.70 


1.2; 












20 kHz 


(fixed) 


(fixed) 












40 kHz 


1.3 : 


• 5:7 '.r 












80 kHz 


(fixed) 


(fixed) 


10 


- 3. 


5 


Koi 


ciV 


0.2 


2.6 


2.4 








Kio 


c 2 V 


0.2 


1.2 


1.2 








K02 


c 3 V 


: 0.2 


1.6 


1.9 








K n 




40 kHz 


1.4 


1.8 








K20 




. 20 kHz 


(fixed) 


(fixed) 












60 kHz 


(fixed) 


(fixed) 












10 kHz 


3.6 . 


• 9.0 












60 kHz 


(fixed) 


(fixed) 








t 


. qB3 


70 kHz 


(fixed) 


* (fixed) 


10 . 


2 


6 


1 

Koi 


c{V 


0.5 


1.1 . 










Kto 


C2V 


0,5 


1.0 


5.0 • 










* 

<1ai 


60 kHz 


0.70 


- 1.4 • 








Kn 


qA2 


20 kHz 


1.3 


. 4.5, . 












40 kHz 


0.90 


2.9 








K21 


qB2 


80 kHz 


0.55 . 


1.5 



i 
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Eiror values are calculated from the scattered results of analysis applied to a series of 
simulated data. 





In Table 4 the relative deviation of mean values (i.e, bias) of estimated parameters are 
presented for 2D-FDDA (present invention) and 2D-MAFID (prior art). In each case bias is 
determined from analysis of . a series of thirty simulated random distributions of count 
numbers. Three cases were analysed. In the first case, models used in data simulations and 
data analysis were identical. In the second case, the distributions of count numbers were 
simulated assuming that particles of the second species are not equivalent but being 
distributed by their individual brightness with a relative half-width of 20 percent. This 
phenomenon was intentionally ignored in analysis, however. Of course, applying a slightly 
inadequate model for analysis produces bias of estimated parameters. The third case is similar 
to the second one except the relative half-width of the individual brightness distribution of the 
second species is 50 percent, which is a usual value for vesicular preparations. It is worth 
noting that methodological deviations are noticeable when mapping weighted residuals of 2D- 
FIDA in cases two and three, but 2D-FIDA still returns meaningful results. 2D-MAFID is a 
more sensitive method to model roughness. 




'8 
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Table 4 



Parameter 

.r * j 

specification (selected) 



True value Bias of2D-FIDA, percent 



Bias of 2D-MAFID, percent 



Case I Case 2 Case 3 Case 1 



Case 2 Case 3 




ciV 



i.O 



0.05 



20 kHz 



+0.1 ±0,2 -0.35 



+0.6 ±0.4 -2;0 
-0.1 ±0.2 +0.8 



-10.6 



+1.8 ±0.6 -15.6 -69 



+5.0 -0.3 ±0.3 +10.0 ; +58 



100 kHz - -0.3 ±0.4 -0,7 



1 kHz 



-1.2 ±0-5 +2.6 



-11.2 



+20.6 



-019 ±0.3 +5.5 



-6.7 ±3.0 +75 



+41 



563 



200 kHz -0.3+0.4 +1.2 



■1.3 



-0:9 ±03 +11.6 . +99 




BRIEF DESCRIPTION OF FIGURES 

Figure 1 shows a count number distribution obtained for a solution of the dye 
tetramethylrhodamine (a) and residual curves corresponding to different fitting procedures 

. * J ' . " 

Figure 2 shows simulated distributions of the number of counts for a case which models a 
binding reaction of a labeled ligand to vesicles . 

- ■ " " t * , .. * '4 

\ ... . * > t t ' ^ " h ' . < * 

Figure. 3 illustrates theoretical errors of the estimated parameters e and ? of a solution of 
single species, depending on the value ofq. ^. 

•* * ' * 

Figure' 4 illustrates theoretical errors of the estimated parameters c and q of a solution of 
single species, depending on me value of c: 

' ' ' ' ' ■ - '■ r / , - * '■ , 

1 * ' - ' 

1 ■ ■ , ■ t * - . " L '" 4. 

* ■ » ( * J 

Figure 5 illustrates theoretical errors of the estimated parameters of c and q of a solution of 
single species, depending on the value of T. 5 



i 
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Figure ~6 illustrates theoretical errors of the estimated parameters c and q of a mixture of two 
species, depending on the value of T. 

Figure 7 illustrates theoretical errors of the estimated parameters c and q of a mixture of two 
species, depending on the ratio ofg^toq^ 

Figure 8 illustrates theoretical errors of the estimated parameters c and g of a mixture of two 
specie^ depending on concentrations. 

. . .'. • . . ■ 

Figure 9 illustrates an apparatus adapted for use in performing the method according to the 

present invention. 

In Figure 10A, the calculated photon count number distributions, P(ri), for cases having 
identical mean count number n but differing by the composition of the sample are plotted. 

Figure 10 B illustrates the distributions of the number of photon counts of pure solutions of 
rhodamine 60 (Rh6G) and tetramethylrhodamine (TMR) ? as well as a mixture of these two 
dyes. 

The exemplary residuals for Rh6G are shown in Figure IOC. 




Figure 10 D shows the results of an ITR analysis (inverse transformation with the help of 
linear regularization and constraining concentrations to non-negative values) applied to the 
curves of Figure 10 B. 



Figure 11 shows an ITR analysis of hybridized (A-C) and restriction enzyme cleaved (D, E) 
labeled oligonucleotides. The curves result from a set of 20 individual 10 s measurements 
which show variations among each other of less than 10%. 



Figure 12 illustrates hybridization and restriction enzyme cleavage of different combinations 
of labeled and unlabeled oligonucleotides. 



Figure 13 is a graphical 
measured for a solution of TAMRA. 

• ■ _ ; ■ ■ < , ' ' " ...■* ■ - ' ' • . ■ ^ " " . \ . ■ 

Figure 14 is a ^aphical presentation of weighted residuals of a joint distribution of count 

numbers, obtained from a mixture of TAMRA (5%-^ 

RRX (rhodamine red X). . v '\ _ • '; .; • .. 

" / - » . ' *r , ' _ ' . - '• * ,.- ■ ' . ' ' ' . , ' . ■ - i ; ■ ' 

Figure 15 illustrates joint distributions of the numbers of, phbtoncounts for the "parallel" and 
-perpendicular' polarization components of fluorescence measured for equal theophyllim 

concentration 

Fisure 16 illustrates the results of 2D-FEDA applied to data from samples of equal 

— antibody concentration. 

s". . • ' ' . .... * ' , '. ' -'- . • - 

Figure 17 illustrates the principle of a two-color 2D-FIDA "experiment with multiple binding 



Figure 18 illustrates joint distributions of the numbers of photon counts in "green" and "red- 
measured in conditions of high (A) and low (B) degree of binding of SMS-14 (somatostatin- 
14) to SSTR-2 (human type-: 




Fisure 19 shows the competition cui>e of uie bmdirig reaction of SMS-14 to SSTR-2 . 

9 : . ■ ■■ r 

?. > ■ ; . . • - :^ • • . ... . _ .• • • ' ■ 1 • ; ■ ■ .. : • 

Figure 20 Illustrates an embodiment in which numbers of photon counts M subject to 
determination of a distribution function ?W in step b) are derived from numbers of photon 
counts in primary time intervals {Nj} by addition of numbers of photon counts from primary 
time intervals according to a predetermined rule. 

Figure 21 Ulostrates a further embodiment in which numbers of photon counts {*} are 
derived from predetermined primary time intends according to a rule in which primary time 
intervals are separated by a time delay. 
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Figure 22 shows a schematic diagram of binding of 3 fluoresoently-labelied molecule to 
receptor-bearing vesicles as observed utilizing the present invention. 

Figure 23 illustrates the ligand binding characteristics of TMR-EGF to membranes prepared 
from A431 cells. 




Figure 24 shows the binding ofTMR-EGF to membranes prepared from HCT cells.. 

Figure 25 illustrates the pharmacological profile of the EGF receptor. 

■ ■ . . . • ■ \ 

Figure 26 illustrates the ligand binding characteristics of Bodipy FL-CGP 121 77. 

Figure 27 shows the inhibition of Bodipy FL^CGP 12177 binding in 40 yl and 1 jil assay 
volumes. 




DETAILED DESCRIPTION OF FIGURES 

Reference is now made to Figure 1 which shows a count number distribution obtained for a 
solution of the dye tetramethylrhodamine (a) and three residuals curves corresponding to the 
best fit obtained with Eqs (36), (37), and (38) below (b). The concentration of the aqueous 
tetramethylrhodamine solution was about 10" 9 M. Primary time intervals were 40 The data 
collection time was 60 s. 



The most widely used spatial profile model in FCS is the three-dimensional Gaussian profile 
with a single parameter of shape, the axial dimension ratio in longitudinal and radial 
directions. Residuals curve with open circles illustrates the fit quality obtainable with the 
Gaussian profile. There are large systematic deviations in residuals. What is a sufficiently 
flexible model for fitting FCS data has turned out to be a rather inflexible and inadequate 
model for FIDA. 



A model of the sample profile which has yielded a better fit of the measured distribution 
^(n) is Gaussian-Lorentzian (open squares), but this model still lacks flexibility. According 



to Eq. (8), a certain function of the spatial brightness B is integrated over the volume; In other 
words, it is a relationship between B and V, characterizing a given spatial brightness profile in 

.1 t * - 

FIDA. For example, the Gaus sian profile yields the relationship 



dV : r ■ ■ ■ - v. . . : / ■ . - (36) 

ax ■■ 



where x = - in B. The Gauss-Lorentzian profile yields the relationship 



:M Jsinh|r. 



^ L-u£ (37) 



Both of the relationships are rather inflexible, i.e., they do/not provide any spatial shape 
parameters to adjust the theoretically calculated distribution to fit the measured data. 

When looking for sufficiently flexible models to fit experimental data, it might be useful to 
apply the following relationship: 

dV ■ . \ 3. (38) 

tit .. . . - .. . • . . • , • - ., • "' ' * ■ 

* . ^ ■ 1 i .... . y 

There is a formal rather than a physical model behind Eq. (38). The fit quality obtainable with 
Eq. (38) is illustrated by the filled squares curve. 

i- - i 

Referring now to Figure 2, simulated distributions of the number of counts are shown for a 
case which models a binding reaction of a labeled ligand to vesicles. Primary time intervals 
were 40 ^/ values of the spatial parameters of Eq. (38) a, = -0.4; a 2 = 0.08, background count 
rate 6 = 1,0 kHz, data collection time 4 s. Curve Jigand" corresponds to a species of c = 6.0; 
q = 6.0 kHz/particle; cx, - 0. Curve „vesicies" corresponds to a species of c .- 0.05; ? = 300.0 
kHz/particle: cr<, -150.0. Curve ^mixture" corresponds to their mixture. Concentrations and 
specific brightness values have been selected to model a characteristic situation in drug 
screening. Fitting of curve (c) returns the values of the five parameters characterizing the 
given ..sample" with statistical errors which are mostly between .3.5 and 6 percent, except the 
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error of er q of vesicles which is 13 percent If. however, o q of vesicles is fixed in fitting, all the 
statistical errors are below 4 percent 

For the fastest data simulation algorithm, one may calculate the expected distribution and 
generate a random Poisson number of events for each value of n independently. Generation 
of a random Poisson number is the following: For a given expected value of events E, a 
simulated number of events S is determined from a routinely generated number R between 0.0 
and 1.0 through inequations 

s-i s 

E Poisson (i; E) < R < 2 Poisson (i; E) 

As a cosmetic error, the total number of events may slightly deviate from the pre- 

n 

given number M A slower but a straightforward data simulation algorithm is the generation of 
a random configuration of particles in volume elements contributing to fluorescence, the 
calculation of the classical light intensi ty corresponding to the given configuration of particles, 
and the generation of a random Poisson number corresponding to this intensity, as a simulated 
number of photon counts. This procedure is repeated M times to obtain a simulated count 
number distribution. 

Reference is now made to Figure 3 which illustrates the theoretical errors of the estimated 
parameters c and q of a solution of single species, depending on the value of q. The following 
values of experimental parameters were selected: c - 1.0; T- 20 jis; a x = -0.4; a 2 555 0.08; b = 

+ r L ' 

1,0 kHz; data collection time 2 s. 

Figure 4 illustrates the theoretical errors of the estimated parameters c and q of a solution of 
single species, depending on the value of c. The following values of experimental parameters 
were selected: g = 60 kHz'particle; T = 20 jis; a : = -0.4; a 2 = 0.08; b = 1.0 kHz; data 
collection time 2 s. 

. * ■ ■ 

Referring now to Figure 5, theoretical errors of the estimated parameters c and q of a solution 
of single species are shown, depending on the value of T. The following values of 
experimental parameters were selected; c = 1 .0; q .» 30,0 kHz/particle (upper graphs); q = 60.0 
kHzfparticle (lower graphs); a t = -0.4; a 2 = 0.08; b = 1.0 kHz; data collection time 2 s, 



Figure 6 illustrates the theoretical errors of the estimated parameters c and q of a mixture of 
two species, depending on the value of. T. The following values of experimental parameters 
were selected: d = 0.1; * - 2:0; gi = 200:0 kHz/particle; ? 2 " 10.0 kHz/particle; *i - -0.4; <n 
= 0.08; b - 1.0 kHz; data collection time 10 s. Note that the optimal value of T for the 
determination of the parameters of the brighter species is lower thanthat of the darker species. 

Reference is now made to Figure 7 which illustrates the theoretical errors of the estimated 
parameters c and q of a mixture of two species, depending on the ratio of q ? to q,. The 
following values of experimental parameters were selected: ^ = 50.0 kHz/particle; c t = c 2 - 
q .5. ax « -OA; az = 0.08; b » 1 .0 kHz; data collection time 40 s. 

+ 1 ' ' > ** * i ' ' " - * 

Figure 8 illustrates the theoretical errors of the estimated parameters c and q of a mixture of 
two species, depending on concentrations. The concenlrations were changed synchronously, c, 
* ti. The following values of experimental parameters were selected; ?i =75.0 kHz/particle; 
<fe = 25.0 kH2/particle; a t = -OA; * 2 ~ 0.08, b = 1.0 kHz; data collection time 60 s. Note that 
an optimal concentration exists at about one particle per sample volume. This is generally 
true, except if less than three parameters are to be determined. 

• ' ' S . • . ' - 

Reference is now made to Figure. 9 which shows one embodiment of an apparatus adapted for 
use in performing the method according to the present invention. Apparatus 10 comprises a 
laser 12 which serves as a light source for wuminating the sample by a bundle of coherent, 
monochromatic excitation radiation 14. Excitation radiation 14 is paralleled by a lens 16 and 
reaches a dichroic mirror 20. Preferably, the angle between the optical axes 18 and the 
dichroic mirror 20 is 45°. The dichroic minor 20 reflects the excitation radiation 14 in 
direction of an objective lens 22 having its focus 24 within a sample volume 26. Sample 
volume 26 and objective lens 22 are preferably separated from each other by a transparent 
cover glass 28, e. g by the bottom of a ^ commercially available iruoro-tit^ 
the sample. The sample preferably includes fluorescently labelled molecules or other particles. 
Due to excitation by an appropriate excitation radiation 14, the molecules or other particles 
present in mesample emit radiation 30. Emission radiation 30 passes the objective lew 
reaches the dichroic mirror 20 which is transparent for emission radiation 30. Thereafter, 
emission radiation passes a filter 34 and a collimator lens 36 on the optical axes 32. A pinhole 
38 is situated in the focus of collimator lens .36. Emission radiation 30 passing the pinhole 38 
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reaches a further lens 40 and, thereafter, is detected by the photo-detector 42. Within the 
pathway of emission radiation 30, in particular between dichroic mirror 20 and photo-detector 
42, an opaque means 44 is provided through which a central part of the emission radiation 30 
cannot pass. This central part of the emission radiation 30 stems from areas on the optical 
axes 32 in front of or behind the focus 24 of the excitation Tadiation 14. Only emission 
radiation 30 that stems from the focus 24 or its direct neighbourhood passes the pinhole 38 
and reaches photo-detector 42. Instead of placing an opaque means 44 within the pathway of 
emission radiation 30, the pathway of excitation radiation 14 is also suitable for positioning an 
opaque means 44. In particular, an opaque means 44 can be positioned between laser 12 and 
dichroic mirror 20. Use of an opaque means 44 as described in detail herein improves the 
signal-to-noise ratio. 

In Figure 10A, the calculated photon count number distributions, P(n), for five cases of equal 
mean count rate, n = 1.0 are shown. The open symbols correspond to solutions of single 
species, but with different values of the mean count number per particle qT. The solid line is 
calculated for a mixture of two species, one with qT = 05 and the other with q T = 8.0 . 
Obviously, the curves differ from each other considerably. The important point is that FIDA 
according to the present invention can unambiguously separate the contributions of the 
individual species. 



Reference is now made to Figure 10 B which illustrates the distributions of the number of 
photon counts of pure solutions of two different dyes, 0.5 nM rhodamine 6G (Rh6G) and 1.5 
nM tetramethylrhodamine (TMR), as well as a mixture of the two (0.8 nM TMR, 0.1 nM 
Rh6GX The main equipment is a confocal microscope (ConfbCor®; EVOTEC BioSystems 
and Carl Zeiss, Germany) routinely used for fluorescence correlation studies. An attenuated 
(to about 800 fiW) beam from an argon ion laser, wavelength 514.5 nm, is focussed to a spot 
of approximately 0.5 jim radius, which is twice the size of a spot in usual FCS experiments, 
and results in a diffusion time of approximately 200 jxs for rhodamine 6G. The excitation 
intensity has generally been kept lower than or equal to a level characterized by about 15 
percent amplitude of (he triplet term of the auto-correlation function. Fluorescence emission is 
detected through a pinhole on the focal plane of the microscope using an avalanche 
photodiode detector SPCM-AQ 131 (EG&G). The distributions of the number of photon 
counts were measured at T = 40 (is dwell time. Data collection time was 50 s. These 
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■ U -rf^ a ^fnr FIDA The results of a multi-component fit analysis are 
distributions serve as input data for l-UJA. ine rcsu. 

given in Table 5 below with * 2 values calculated according: 



(39) 



^ n/te to iWh of,1hi di^Won of fit 

In said niul«nen< analyst, on. » *. ,M dtsW ^on of 
„*. of photon counts, L**, a certain nuntb* of fluo^cen, spec,.,, and estnnaes 
a.too^n'wncOTlrauon and spwiflc brightness values. .■: •;■ ; ■ : • 
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TableS 



Sample 




A3 


C 


g(kHz) 


x 1 


Rh6G 


-0.3 8Q ±0.009 


0.077x0.003 


0.461 ±0.003 


107.2 ±0.* 


0.97 


TMR 


-0.427 ±0.032 


0.08410.014 


1.517 ±0.012 


36.56 ±0.29 


0.8] 


Mixture 


-0,380 (fixed) 

— . : , 


0.077 (fixed) 


0.103 ±0.012 
0.738 ±0.011 


109.1±4.0 

* * 

37.4 ±L0 


0.77 



In all cases, the background count rate was fixed to the value of 1.05 kHz, as measured with 
de-ionized water. a% and a% are pre-normalization values when 0\ is fixed to 1,0/ The 
theoretical statistical errors in Table 5 correspond to theoretical weights: 



f- (40) 

This formula is derived under a simple assumption of ,V independent measurements of the 
number of photon counts. In reality, consecutive measurements are correlated, therefore the 
errors of estimated parameters returned by the fitting algorithm underestimate real statistical 
errors. We have empirically determined from a separate series of 30 to 200 measurements that 
these statistical errors are greater than theoretical ones by a factor of about three. 

The exemplary residuals for Rh6G are shown in Figure 10 C. 

The results of ITR analysis are shown in Figure 10 D expressing the distribution of the 
number of particles as a function of their specific brightness. Such ITR analysis is an inverse 
transformation realized with the help of linear regularization and constraining concentrations 
to non-negative values. ITR is a valuable tool especially for samples from which either no a 
priori information about the sample composition is given or where the sample composition is 
heterogeneous. The dashed lines correspond to the solutions of the single dyes (Rh6G and 
TMR) and the solid line to their mixture. The ordinates give the mean number of particles 
within the confocal volume element (left: single dyes; right: mixture). For equipment details 
see description of Figure 10B. The ideal outcome would be single 5 -peaks for the solutions of 



.A 



single species and two S -peaks for the mixture. In reality, the width of the ITR output spectral 
peaks is determined not only by the true width of the distribution of specific . brightness values 
but also by the accuracy of input data and the particular realization of the Unearregularization : 
simply, if a broad spectrum fits experimental data as well as a narrow one then ITR prefers the 
broad one. ■■ ■ ' '.. ; , •V ; '\"' ', 

As a further example, the method according to the present invention has been applied to study 
the hybridization of 5K6^boxyterramethylrhodamiiie (TANmA)>Iabeled 40mers with 
either labeled or non-labeled complementary oligonucleotides and me subsequent symmetrical 
cleavage of the DN A hybrid by the restriction endonucleases Hind III and KpnJ. The specific 

oligonucleotides used in this study were - 
TAMRA-AAGAAGGCKjTACCTTTGGATAAAAGAGAAGCTTT^ 

OIigoA)and 

TAMRA-ACGGGAAAAGCTTCT 

OligoB). . ■'. \ . " .. 

They were purchased in HPLC pure quality from Applied Biosystems (Weiterstadt, Germany). 

All measurements were carried but on the above described FCS reader at excitation/emission 
wavelengths of 543/580 nm using a 4 mW helium/neon laser (Uniphase) attenuated to 
approximately 300 uW. The water background was found to be below 800 Hz. For the 
measurements, sample aliquot? were diluted to I nM and 20 ul were assayed in a 8-well 
chambered coverglass (Nalge Nunc) at room temperature. The hybridization reaction was 
performed in 70% fbnnamide contairung lOmM Tris/HCl buffer (pH 8.0), ImM EDTA, 
0.2mM NaClandan oligonucleotide concentration of0.5uM. Denaturation was at 95°C for 2 
min and subsequent hybridization was at 55-60'C for 40 min in accordance with the 
optimized temperature Tx which is 10-15 degrees below the melting point Tm (Heating block/ 
Techne). Tm was calculated as follows: Tm - 81.5 +16.6 (lo g 10[Na+]) + 0,41 (%G+G) - 
600/N (N = 40; %G+C = 42.5). Restriction digest analyses of the hybrid DNA were performed 
by the restriction enzymes Hind III and Kpnl. The restriction site was chosen in order to 
obtain fragments of different size. The cleavage reactions were performed in 3.3mM 
tris/acetate (pH 7.9), 1 mM magnesium acetate, 6.6. mM potassium acetate and p:imgVml 
bovine serum albumin at 37°C for lh. The reaction course is characterized by significant 
shifts of the fluorescence intensity per molecule within a range of one order of magnitude 
(Figures 11 A - E). It is important to note that the resulting double product peak corresponds 
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to one single hybridization species in which no none-hybridized species are detectable. This 
could be easily demonstrated in hybridization/cleavage experiments using labeled/nonlabeled 
oligonucleotide combinations where signals with different brightness compared with the 
double labeled hybrid have been obtained (Fig. 12 A3)- One explanation for such a double- 
peak structure could be conformational fluctuations resulting from transitions between two 
conformational states. This behavior is described using a three-state model of the 
conformational dynamics with a polar, a nonpolar, and a quenching environment of the label 
(Eggeling et at, Proc. Natl. Acad. Sci. USA 95: 15S6 - 1561). From the present data it can be 
concluded that the conformational fluctuations are significantly dependent on the DNA length 
and the number of labels. The small digested DNA fragments exhibit a single brightness and 
the concentration ratio of the educt double peaks shifts toward one component when one-label 
hybrids are studied. As a further result, small fragments show low molecular brightness and 
large fragments correspond to high brightness values (Q2 « Qs, Qi « Q4)- 

Based on symmetrical cleavage sites for both enzymes DNA fragments of equal size and 
subsequently very similar molecular brightness are to be expected Figures 12 A and B 
clearly illustrate that such a result has been observed and the theoretical relationship Qi + Q 2 = 
Q3 + Q4 is experimentally confirmed. Restriction digest analyses of the hybrid DNA were 
performed with the restriction enzymes Hind ID and Kpn I either alone or in combination. The 
restriction sites lead to symmetrically cleaved fragments (2x12 mers, 1x16 mer for double 
digest). [•] corresponds to labeled DNA. Qouo^ intensities of DNA hybrids, Qj: intensities of 
cleavage products. The bulk of all the measurements of cleavage products using single and 
doubly labeled DNA structures makes it possible to classify all observed intensities with 
individual structures in a single measurement - a property of the method of the present 
invention which is not possible using any other known analytical method. 

Reference is now made to Figures 13 and 14. As the central optical part of the equipment for 
2D-F1DA, a confocal microscope is used, like in FCS spectrometers (Koppel et ai., Biophys. 
J, 16: 1315 - 1329, 1976). For excitation of fluorescence, a team from an Ar or green He-Ne 
laser is attenuated by neutral filters, passes a beam expander and is directed to the microscope 
objective by a dichroic minor. In a number of experiments with slowly diffusing particles, 
beam scanning in combination with sample scanning is used, as a tool known from laser 



scanning microscopy. Fluorescence is collected % the same objective through the dichroio 
mirror, and is focussed to a corifocal pinhole which serves to reject the oiit-of-focus light. 
Resolution in the longitudinal direction is additionally improved by using a concentrical 
opaque Spo t dosing about a quarter of the aperture in the fluorescence collection path of the 
microscope (see description of Figure 9). The light which passes the pinhole is divided by a 

• ' * > 

beamsplitter for detection by two detectors. Depending on the general type of a 2D-FIDA 
experiment, the beamsplitter is either a polarization cube, or a dichroic mirror. In the first 
case, a common spectral band-pass filter is used, while in the case of two-colour FID A, each 
detector has a different band-pass fitter in front of it The photon counting detectors are silicon 
avalanche photodiode modules SPCM-AQ-131, EG&G Optoelectronics, Canada. The TTL 
pulses from the detectors are counted by a two-channel counter, constructed by EVOTEC 
BioSystems AG (Hamburg, Germany) as a plug-in card of a computer. The count number 
distributions are calculated during reading data from the 32 MB onboard buffer, which itself is 
an online process. By feeding the detector outputs to a correlator, FCS measurements can be 
performed in parallel with FID A experiments. 

. . _ ^ ■ i J ' . ' ' 

*■ » ' . » • . i < ~ » _ . • - . ■ 

The levels of background count rate for both detectors are determined by a separate 
experiment on bi-distilled water: The main contributor to the non-fluctuating background light 
intensity is Raman scattering from water. 



The radius of the monitored sample volume can be adjusted by selecting an appropriate 
expansion factor of the original laser beam. The focal beam radius of about 0.6 um is used 

-j 

yielding diffusion times for simple organic dye molecules (e.g., TAMRA) of about 260 us, 

In. 

which is considered long compared to the 40 us dwell time of counters, so that the 
approximation of constant molecular brightness Airing the counting interval is valid The 
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excitation intensity is adjusted as a compromise between a high count rate per molecule and 
low population of the triplet state. The triplet , state population was kept at about 15. percent 
level. Higher triplet population values might significantly distort the apparent, spatiaJ 
brightness profile. The population of the triplet state was measured with the correlator (ALV- 
5000, ALV, Langen, Germany) during the set-up of a series of experiments. 

For methodological test experiments, two different dyes were selected, carboxytetramethyl- 
rhodamine (TAMRA) and rhodamine red X (RRX). These dyes have different emission 
spectra, as well as different extinction coefficients at the laser wavelength of 543.5 nm. In 
these experiments/ a wideband 40/60 beamsplitter was used in front of the detectors. The 
spectral filter of the „red 6C channel has the center wavelength of 605 nm and F WHM of 50 nm 
while the corresponding figures for the 3 ,yeilow-green" channel are 575 nm and 30 nm. The 
dyes were diluted in distilled water so that the average number of molecules in the observation 
volume was in the range of 0.5 to 2.0, corresponding to concentrations between 0.23 and 0.92 
nM. 

For each experiment, about 20 \xL of the sample solution was placed on a coverslip separating 
the sample from the water immersion objective (Zeiss C-Apochromat 40x1.2 W Korr). Each 
distribution was collected for 60 seconds. Each of the dyes was measured separately, but also 
mixtures of the dyes with different concentration ratios were measured- The parameters 
describing the spatial brightness profile were determined from an experiment on TAMRA and 
were fixed in subsequent analysis of other samples at values a x ~ -0.405 and a 2 0.0772. 
Mixtures were measured and analyzed in order to see how well the method returns values of 
specific parameters of the two species. 



TAMRA solutioiL Va,ues ceding" the ^ are logarithm values of** 
Wobtain a giVen pair of count numbers and *, Fi*^ of *. disuihunon returns m*n 
n umber of panicie, cK = lAW^^M***** » Kd " ^* " W - 4 

-0,*,, and for the „ y e ilo ^on»n e l „ - 50.9 i<j.3 kHz. In Table 6,^ results of 
^ of test experiment ^ present 8-** have no. been specified by 

oMep^on^ausorpucm of uye.olecu.es. 

„ ^ass Places fences W cor^rauon, significant Because of the same 
eon^on values are no. v~U produced W s^np.e i sample; a sMft of concent 
vaK.es to realization «o realization is a!so sometimes observable. Specific bright values 
are well reproduced, however. Values of statistical errors present are theoretical values 
^spondU* M d»»ssun»»enofuncorre,a«d measurements, mutepHed by dnee, which is 
an empMcal factor. In addition to: aatistical errors, deviation, ftoro sample to sample of a 



modest size are also noticeable. 



Table 6 

arocnber molecules per in „red», inflow-green, 
sample volume, cV , ; 



ch, kHz 



TAMRA 


1 


1.128 ±0.005 




2 


1.139 




3 


; .1.160 




4 


1.171 


RRX 


1 


1.892 +0.009; 




"\ 2 - 


1.859 . 




3 


1.850 



79.6 ±0.4 51.3 ±0.3 

79.4 S0.9 

79.4 50.8 

79.0 50.6 



48^03—^" 15.63 ±0.09 
48.0 15.63 ^ 



1011 



47.3 
47.9 



15-44 
15.58 
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TAMRA and 1 
RRX * 




TAMRA and 1 

RRX, new 
sample 2 



4 



# 



: 23 -.12 -99 ; Mm&m^mm%P 221 134297- 
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0.9910.05 


78.0 ±1.0 


< 1 O -Li 1 


1 ■ 14 ±0.05 


49.1 ±1.1 




1.01 


76.8 


* ? 1 


1.08 


49.5 


15 1 


. 0,99 


77.9 


J 1 .0 


1.03 


49.3 


1< 7 

O.i 


0.97 


76.9 




.1.03. 


50.0 


O.J 


0.93 ±0,07 


74.6+1 3 

* " J ,* 


49.0 ±1.5 


2.51 ±0.07 


476+0 8 


14.7 ±0.5 


0.98 


75 2 


49.2 


2.38 


46.8 


142 


0.96 


75.0 ,. \ ,\7 


49.4 


2,25 


48.0 


1 : 

14.7 


0.94 


76.1 


49.9 


2.16 


47.6 


14.8 



to. « p™ ^ ns „„ ,4 vi sualiM3 ^ ^ of ^ a 
™ wi<h a of TAMRA ^ RRX m ^ ^ „ ^ ^ 

-fr* when two Spe cia, ^ ismci .„ ^ ^ ^ ^ ^ 

Fe.™; «fa_. b=m«„ ... mjmMd ^ ^ 

i» M, ca*. , h ooum mate <*M W ». „ K<P . oott*,, «, w is „ 60lml nilfflber 

obtained by the „yellow-green" detector. 



b Mfad Wscaocc poM^a ^ avarage intensities ' of two potent 
of fl " OTK »°- - ^-MftM. of*, intaosito a« no, of 



fluorescence polarization value, of a s amp* oon^.TO^^^g partner 

'^■^i^J-i^iw^. canaisob^rfom.edta real-time, allowing 
thaiiianalysis it ^ociation and dissociation reactions. One af most widely-nsed 

piLii-iW W^-f^ii .^<^--* 1i * i,,1 .^:■**" a, * i,,^0 ?■ : 

ical immOn^ssays for. «A#* . V decade- 'ti, B**.0*^ 
poison in irnmunoassays) has we»ed Vantages over conventional heKrogeneous 

investigated because .he separation 6i^^ f^^-^^f^-^'- 
ligana binding curves ^a-^*^**- of poiari^on. meaning to. 

^nistic ^i^^*"^-^?*** 

volume as well as in mass restrictions. \ 



With a dimensional fluoricenc* -t^lJipi^.^"^* 

invention, ihe * floOKS ^ 

utiiizal, thereby overcoming the talons mentioned above. 2T>FIDA direcUy determines 
: two sr.«.fie quantities pe, ^toAeence specie^ one measurement: the ft,o^n« 
intosity per n,olocu.e and th. anisotropy of a given model. Based on to supplementary 
information, Ihe delineation of all pamctpaing ^ and ever, the quantification of tie 
bWing behavior * possible. 2D-FIDA anisonopy is an ided tool for & quatnltative 

^on o f ^•*m*^^:^ : *r* i * * nd muW,nKi2 *" 



phenomena- 



~ TM^ te the ^ ^ ^ ^ ^ ^ ^ 

and wthout an incorporated spacer • In classical FP analysis the- . .: 

analysis ttese conjugates exhibited low 
amsotropy values (0.03? and 0 055 iesi>Mii»i„> . • . 

ftey a ciW basis fcr „ 

For binding^pertam^ ^ M]Mto n . 

soTOons of anbbc^y a^ antisens wne diluted in a 

» ^ — - — „ incubafcd 

* 30 mmte . ^ ^ ^ ^ ^ ^ ^ conjugates the 

measured joint distribulicms of photon ™,™ „„_•. , 

Photon count numbers and results of 2D-HDA applied to 

.samples at different antibody diltnion values but a constant ligand concentration of {£} • 2nM 

^^^^^^^^^^ 

««mes a high number of receptors v^,vi M ;„ . 

receptors, ves lc Ies m samples with a low binding degree can be 

d^guishwl fiom vesicles in samples with hi<* w a- , 

^^ *&**dmg degree by z ^^yu^ 

0fV ^ * ^ * ^ « ^ but ^ , ^ fc a 



ta ' W »»- Co— » *on, two f.^. , 

in M «««*.» . 

• w^mP of TAMRA-labeled sbmatostatin-14 

. _ Trt-psent invention* the diucuhb wuv* . 

analysis according to toe present mvwm » 

t . ^ ■> w B h afrMtV somatostatin receptor SSTR-2 (P. 
(SMS14-5TAMRA) to the human type-2 high attmity 

, « T PWii 289 ^ 
Schoeffter et al., Eur. J. Pharm., 28*, leo 

rffcv CCL39h SS t2 cells. Vesicles were ^ pr^d from this ceU luie 
The receptor was expressed by CCL3ynss^.c C1 

HBPBS (pH : 7* 0.0,««*> fluortisurfectent FC- 1 35» ,33,,PMSO in *, 

^ofp^^-A"^ 
in Figure 17. / : 



K^Vas - .» ^ * * - *» *'* ^ " 

^pa^ouBha^^ 

the collected ftiorescencc beam passes a smiric v 

^ on. «*. ~» -C^ 6,0 * BWHM" « - ** 

^ n^ fl~ce *. *. «*o .abeU Becan* v^es 



« ^ scan^of M Hz ^ 100 , _ ^ ^ ^ 

scanning of 50O.(jjn pet 8 s data collection tune io ike wber directions 



a ^ de^ „ f „ coo^. „ ^ „ B, * veacle ^ ^ (SMS]4 . 

>TAMRA) „ competed off* „ ^ of , ^ ^ rf ^ — ^ 
(SRIF-14, SIGMA; , * con^to, 3 nM , is ^ ^ ^ 

fl*^ of te taW .o ^ mt ^ !S , ^ fc ^ ^ ^ 
^^^^^^ ^ to ^ ^ ^ ^ ' 

Ob, . ^ mpMra , fi , ^ ^ ^ — ^ ^ 

and „ Wence W6toMS ^ ^ ^ ^ 

tovesie.es ra ^ ^ yMjclK „ ^^ — ^ ^ ^ ^ . . 

green fluorescence of the vesicles and an incrai* i„ . 

an increase u, the rancentnWcai of See Ug^. m ^ 

ttat the lijand does not completely dissocial 

era. assoc.™ Son, ^ VMlcfe fe ^ . 

Hex, fc eontpetito, has been ^i*,* ^ pM (w ^ 

». a** ana^s ieads » a ^ ^ MBpelidM _ ^ ^ ^ ^ 



in 



values 



hi* and low W^t* -* «" ^ ^ °" 
* t es».ved with 2D-FIDA analyS * * ^ i! iD ^ 

if one aims a, me gallon of assay — Chi*") * <~< 

^ «■ -*> W • ****** KUab,e r-f * 



screening* 



photon courts W subject * Ration of V diaribunon mn*on *0 ■.- ;*» » - 

■ -u <rf by edation of numbers 

derived from numbers of photon counts in pnw»? . 

d L«o„ Action ft.) are : ^ from *e u™»ers of photon «~ '» 

. ^ -Vtv where Mis an integer number 
^ time intervals {Nj> according » the rule n, = . 

easing how many times ft. <i»» inter^ in which \.,> is denned M-tM*£ 
X te^J Line N shows the prima, ^interval window, tine , — me 
prtonrj- time windows chosen to calculate n, according .0 <he rule. v - 



- " .i- ■ , -..1 . r 

21 shows a tather entbodunent in which numbers of photon a» ^ 



from predetennraed primary time u«ervals aocor^ • . .... 

* j^ted by a time delay, . ^cblar, * 
n. -t(*~^**~ M i - Po^ve mteger numbers, **> are numbers 
of phi counts subject to de—on of a distfbuuoo funcnon *n> , and («, are the 
numbers of photon counts in ptimao- "me intervals. 

Reference is now made to figure* 22 to 27. . . . x ^ciplines of 

The study of function and activity of cell surface receptors « central to all discipl 

modem biology and pnarmacowg? ^~ nuta; signals across the plasma 

within a muW^Mar organism since they transduce extracellular signals ac 



- o Prcteta ^ w Keptors (GPCR) ^ ^ ^ ^ £ 

^ a "°"" ^ ** * - *h' ceU cation m 

Biopolymers, 43, 339 - 3«6, 1997: Riese DJ & Stem D V « 

Pereson H. & Ibanez C F Chit o»l n N. „i „r . '"™' 

■w to. Opta Neurol. Neurosurg. 6, II ; - ] S , 1993). O pwem- 

<— »e fa*, ^ ^ tmm ^ m> ^ 

P^GSste by binding extracellular factors such as iw^tranamitMis and I»nn6nes (Ji et J 

1998- WessX., FASEB J. 11,346 - 354,1997); ' • \ 

b most oases, seep** are studied biochemical* and p^logMy using a ^ 
^ or » *,e cCs is monito^, ^ ^ ^ ^ ^™ 

»** m be appbed ■ .^W^ W W ^ te 

fluoreseen: ]i 61n d wjft fluorescence into.*, ^ • ^ 

„. "uuiescence tmcra.iy distnbuoor. analysis rnnaV 

accontag to the ,»eseiit invention, these m^~- , J 

conce.**™ of fee aid WW- oetermmed ta 4 siiS 
mra ^« »«K,utu« need for separettonSmceFtnA,... , ^ 

Wtbe de« M volum, is ^ ^^?T'' a 
VreoWtni i \. i "W W sea ofl femtoliter. The Staple volumes can 

. Ml wbtch c^nUy fe * ^ |imit « ^ 

^ Wd^ ^ ■ % re* 



washing stepsli.e. tn c 

signmcant reduction in reagent oonsumpto , 

v -^^«T,t role in a large number of diseases 
• • • • „ We >i^ribed above play an important roie m 4 i«e 

Both classes of receptors descnoea aoo p v • aT ,d Ali)S and areihus 

* 2 j ^;™riv^ diseases cardiovascular disorders ana au^o ™» 

^ as cancer, neurodegenerative ^ ^ - recep tor tyrosine Iqnase 



imoortant drug targets the EGF receptor (JEGFR/afbp U. auu — - 
^^^^^ 

The following experimental protocol has been used: ; v 

^ , ^' ^dled epidermal growth factor (IMR-EGF), unlabeled 

^ — (Bodipy Ft: M.difluoro-5,7- . 

epidermal growth factor (EGF) and Bodxpy ^ ^ ^ 

chloride) were purauwu ^ mr-n r^lls were obtained from 

<»4111 cells and human colon carcinoma tumor (UCO ° u 

EOF) wre purchased 6»U D S>sttms .^^^ 

^stable clones were selected as V - ^ A M ^j; _ 12290,' 1997). For 

Iate s V., Ilcin U Lahse MI., Proc- M Acad. Sci. U.S.A., 94, 12285 ^ 

• o v l O 8 cells were resuspended m 2 ml ot ice coia 

inhibitor cocktad (Complete Mini, u _ we re removed by 

^ by Bounce homogeni^ 

elation at 900 x 10 min. J***-* hufier and 

% v 70 min A-Q Membrane pellets were resuspended m 0.5 nu binding assay 
.. xg for20mm(4C).M ^J?^ ^ ^ Uemhrax)es 

the protein concentration was determined using 

B ^ ^^.^^^^ Ito afinal voter* of 40 nKsee 
aifereotconcenlranolBoffluCKScentllganoaa r ■. . 20 ^ 

Figure legends for details). EOF receptor (EGFR) binding studies were earned out in 20 
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Hepes P H 74, 140 rnM NaCl, 5 mM KC1. 1.2 mM MgCl 2 , 1.8 mM CaCl 2 , 0.1% BSA and 
0.1% CHAPS. The binding assay buffer for the 02 adrenergic receptor contained 75 mM 
Tns-HCl P H 7.4, 10 mM MgCl 2 , 2 mM EDTA, 0.05% BSA and 0.1% CHAPS 
Measurements were earned out in Swell glass chamber, (Nunc). For experiments ^sample 
volumes of 1 ul, proprietary liquid dispensing unit, and high^ensity .ample carriers 
(Nanocarrier)^ a total capacity of 1.5 (il per well were used (EVOTEC BioSystems AG) 
Data analysis. A confocal microscope was used (Coafocor. EVOTEC BioSystems AG and 
Carl Zeiss Jena, Germany) with an optical set-up described in Kask P., Palo K., Ullmann D & 
Gall K. (Proc. Nat,. Acad. Sci. U.S.A. 96, 13756 - 13761, 1999; the contents of which are 
herein incorporated by reference). The number of photon counts was measured at 40 us dwell 
t,me ^ a data ■#«- of 4 ,6 s per Wjfcwi with TMR-EGF were 

earned out at excitatiotVemission wavelengths of 543/5S0 nm using a helium/neon laser 
attenuated to appro«y 300 uW. Measure 
our at excitauon/emission w^vdengms o^ 

approximately 2 mW. Like PCS, me method according to the present invention relies on the 
«ndom drffusion of fluorescent particles in and out of the confocal volume (Kask P Palo K. 
UilmannD. & Gallic, Proc. NaU. Acadi Sct\ US A., 96, 13756-13761, 1999; RiglerR / 
Biotec^Ml, 177- ^ 1 W 5) : 

compared to single molecules such as proteins (mean, diffusion time t = 10 -100 ms versus 
approximately 200 u, - 1 „ for singIe ^ ^ ^ - - . _VjV ^ . ; 

scanntng approach in order to be able to detect sufficiently huge numbers of fluorescent 
****** vesicles. This is achieved using a laser beatn oscillator and a mobile sample 



stage. 



AT«> I~» of .wo differ, fllmssMnt ^ ^ tos^^ 

**" i- 

pho,on distribution was fit wW, « to ^ p ; ^ 

^ **» mm**?****, comports. <W«« of Wand fi« 
»gaod w*, deduc* ftom tha, dau. Ib. of HgaM bOMd „ v raictes * ^vc <*her m 

ON ligaod bound) or as fhefioaa, bindi^ fc^ ^ ^ , - 
*U t^m wctc calctd^tod from tbe oo^do, ' ^ batf^,) 



binding 





23-12-99 :.- 
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34297-> 
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46 



Y . & Lff W. H. <Bioc„em. PWn^.. 22, 3099 - 310*, ,973). Expenmen*. data were 
analyzed using Prism 3.0 (Graph Pad Prism, Sax Diego, USA). 

F,_e totensi* — andvsis (PrcA)a.,ows tbestodv 

of fluorescent molecu.es if . change bright o«urs upon »** Q-^e 

ai describing mbie^ **** « *— ^ ° f 

fluoresce species in «, sample W ^^respccave absolve concents a« 

fere* shows a schematic diagram of binding of , «^ 
^r^arini vesic.es as observed bv the ,ne«hod according to ft, present mveutton. Tne 
■^^l.'i^W.^^'"**'*" A: In the absence of membranes, 
*. confbcol volnme conmins eiehisively *ee fluorescent Ugand mo.ecu.es. 
brighmess of ■h.'aiU**** an intrinsic . Proper* of tne. fluorophore and me 
»to* ft is a«ached to. B: A brightness diStribu^on obtained for the sampi. to- » 

Ugand binding occurs and the.confccal voiume is popular wim bom veside^ound and free 
„: Brightness di^bnflon.o, sample shown in C. Tte molectuar brtghtness of 
%J .igandW- fl» ^ molecular bright of vesicles w* 98 ^ mm, case 
(,4 times the brighmess of free Ugand, corresponding to 14 ligand mb.ecu.es bound per 
veaele) and depend, on the receptor expression leve, Typi«U y , 10 - 500 hgand moleetues 

.~wto ^ . — * *• «- *—« 

■ does io, change upon binding, bom me molar concentration o, bound as - as me moUr 
conizations of unbound ligand are optative* determined ftom a stogie me^uremen. 
without any separation being necessary. , 

* method acccrtdingto me present invention has been used to eh«a«enzeme mteracnor .of 
TMR-labeied EOF with* receptor, the EGE receptor (EGFR/ErbBl). As ■ 
«Us were used which « W to overcxpress EGFR (Mc * '• 

■ CeH Bio, 30; .265 . .27*,. 1.998). Figure 23 shows the Hgand binding 
TMR-EGF to membranes prepaid ftoin A431 cells. Ligand concentrations from ^ 7 
oM were incubate, with a crude membrane fraction a, a fina. concentration of 50 ug/ml m a 
miai voiume of 40 ui. Afier an incubation period of 20 min at room 
^asurements were carried out, and; ** and bound Ugand concentrauons deternuned, 
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Binding of TMR-IabeJed EGF to the membranes is saturable and can be competed bv 
excess of unlabeled EGF Non-st**;^ k- a- \ "" competed by an 

EOF The i, „ 8 — detenUined in & presence of 1 m 

EOF. The results shown are from a tynical Yf*rim 0 „*. . , . ^ 
■ ma ^n^ents, the standard deviation (SD) is shown 

35 error bar * de ^ed dissociation constant (K D ) of 3 1 + 0 2 a ju - 

easier reports where fl uoresce 0tIy ,abe U ed EGF l ^^T^^ 

^ — <<**W« - * *L Chen, to^TZZZ T 
(Bmax) of 60 pmoi/m. was observed which ' * leveI 

P 0 was observed which corresponds to approximately 6x 10 s bindine site, 

-^^^ , f lt , Bioc&em. Cell Biol. 30, , 2o5 , 2: ' 

P " " * * 

In order to demonstrate that the method accoixlingtomenreseht inv^ 

^ mudl ^ „ pw„,„^.^ ion ^ 

SIT "VT" * of a ° f *■*»■ 

«*« ^ ^ (EOre/ErbB-,, BoB-2, E*B3 «1 & b!W Ri 4se a . B 

chancttri^ -n . „ w will HCT manbfanes has been further 

•toaacnzvl. The n^fe, of maamaijaa polyfmifc tartm ....... . 

6 v* lMB EOF). K 1(W e 25 Xmm, fc f harraaoobg,-^ profj, of a* 



The affmnyof uolabeUed EOF, be*** and **ltf^>*»i* 

EGF receptor. The animiy ui ^ iriru hatinc HCT membranes 

.v, w fHB-EGF) was determined by competition assays mcubatmgm, ^ 
growth factor (HB-fc^r; « ^ ™,»,itnr Binding was determined by 

* tur F.GF and the indicated concentrations of competitor. Binding 
Wth TMR " BG F " ff to the present invention and ^ to ^t mfc, absence of 
^^«^.^^^^, iVb .^ for EGF, 0.23 ± 0.02 nM.for 
competitor. The affinities obtained were 0..5 ±.0.12 nM ... ; , 

betacellulin and 1.4 ± 0.3 nM fox HB-EGF. 

the 62 adrenergic receptorwai » weU laiown model 

system tor u» 6 . ■ , „ 2 anergic receptor antagonist, 

,,<,., .11 1992). Ftoorescentty-labeled Cvir Ul//. » ' " » : 
1587 - 1 1", . .„i _ 0134 19941. Figure 26 

as a Ugand (Heithier et A, Biochemistry 33, 9126 - 913". '» . . 

chafes of Bodipy H*GP UX77. Increase 

determined in the presence wj. r r r , . ^ »m diown as error 

- 1 Q ^„„^mf»nt* standard deviations (SD) are snowni* « 

Bmdms of the Hga» ^ ^ ^ a^daM filter Hading assays 

mritWer etal., Biochemistry 33, 91 26 - y . • ■ ■ 

Sll unfcbeM COP 12.77 ^ a Ki * 0, nM. indica** « «. 

1 ' * " , ' 

u - v.-^n of BodinvFL^GP 12177 bindmg in 40 ^1 
Figure 27 illustrates me inhibition of Bodipy v^w 

F,gure m 77AVwas used as a competitor at the concentrations 

volumes. UnlabeUed hgand (CGP 12177A) v«s use _ ^ 

shown, Bm*ng was no^^ - rflHl 

f^n 1 and 1 til For experiments involving sample volumes of 1 Mi, 
ouchisa.nplevolumesof^ulandl ^^^^^^^^ 

proprietary liquid dispensing units and sample earners were us v. 

The affinities obtained for CGP12177 were 0.7 * 0.1 nM in 40 ul and 0.9 * 0.1 in 1 ul assay 

T^sexperWnten^ , 
volumes. This expenrnent emp volumesandmus 0 f the required biological 

opportunity for a significant reduction of sample voium * 
material for high-throughput drug screening. 
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Studying ligand-receptor interactions .is an tmportant tool for the biological and 
pharmacological characterisation of membrane receptor function. Particularly, G protein 
coupled receptors are of high interest for pharmaceutical drug discovery programs siace 
approximately 30 - 50% of all drugs on the market today are estimated to act on GPCR 
(Gudermann et al. J. MoL Med. 73, 51 - 63, 1995). Traditionally, the majority of ligand- 
receptor studies are carried out using radioactively-labeled ligands, where bound fiorn 
unbound ligand is separated by either filtration or o*tf»*i fa recent years, scintillation 
proximity assays have found widespread use for receptor-ligand binding studies Here, 
bmding events are detected by the proximity of the radioligand to membranes immobilized on 
a vmm-comnrns bead, thus eliminating me need for separation (Hart Mol. Lmnunol 
16, 265 - 267, 1 979). However, the disadvantages of these methods are potentially hazardous' 
exposure to radioactive compounds, limited shelf half-life of labeled ligands and special 
requirements for handling ^^it^j^ a ^ ofQm ^ 
based detection methods have bee* applied to receptor-ligand studies in recent years which 
are based on the detection of total fluorescent (Heithier et al, Biochemistry 33, 9126 - 9134 
1994; inglese et al. Biochemistry 37 2372 - 2??7 lo«wn ni. 

- ,8^37, ,5 8 50- ,s 86 4, l^ or^^e^o^^ ^ 
al. Proc. Natl. Acad. Sci. USA, 11122- 1 127, 1999). 

• * 

Based on th, mefcod „ ^ ^ ^ ^ 

analysu « * » d^*e ouanftatf^ 

teceptots. Since ITOA allows interactions „ be ^ ^ stagls ^ ^ ; 
offers significant advantages over the methods described above. .Using EGF and P2 adrenergic 
receptors as model systems it has been shovrii in the examples above that affinities and 

cxp^ob level, „ be stated h a sitnpfc tnanner. me resnta, 
the rWoIogicd profiles obfcined „sing mA ,„ ^ ^ 

radioligand binding assays. It has also been shown that FFDA can be used to study tigaod 
m*~ interns aaog cells express** te w leveis of receptor. This is of parties 
for OPCRs since ^ are typical* fonnd in low density in prinury ^ ^ 
overexptession of facta! receptors - a ^ ^t „ ^ ^ ^ 

the increase in panMe bright ^ billdillg of mow,* auo^ to 



vesicles no changes of fluorescence parameters s ^ as fluor«cer,e ,^chir*« spectra! 

z% — » °- ^ ***** **- . 

colons ti^i*.***^--*^*^ •"7"^^ ' 
providing an Wenud *on«ol <° -* ""^ 

significant advantage of HDA is *. ^ vdumc of detection, w hich is approxuna^the 

si of a bacterial ceU. Minianrr^o. of assay — is . 
ph^aceutica. drug screening since, tocreastag numbers of chemics! compounds have to>e 
Lreened agair^ a growing number of bioiogioa. targets, caption of brrtogrcal ^ 
3^ of chenuea. or neural i^^^i^^™*** 
preferably confbc^rhods based on the present invention are *erefor* rdeally 
Tese purposes since the sample volume can ,» .educed to i rJ._-W._i. current hm« of 
Unhanding, without compos the signal-to^ ratio, therefore offenng «u,ue 

TOA to been succesafidly appied to study a number of diff«en. W *s of O protern-co-pled 
receptors in our laboratory. These include rccepttts, o^arnir, receptor - mt a 

number Of pepnde r^. T^ ^ ^*^ ** «* * ^ 
moleiufcs (agonists or antagonist such as described in this study), and pxc^ns (data 

n0 , shown), to addition, the technology can be successfully applied to measure butchng of 
molecules to li«e cells in a, U^T*^*^'*'^?^?^. 

tissues directly in high-throughput screening. • 
la conclusion, FIDA is a* idea! method for ligand-receptot studies: It offers .ginfican 
Stages for research, high-throughput screening and for pharmacological profiling since it 
■d.o^.ia « of radioactive isotopes and of separation steps, isy^*^*** 
allows a significant reduction of reagent consumption. 
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CLAIMS 

t 

1. A method for characterizing fluorescent molecules or other particles in samples comprising 
the steps of: 




a) monitoring fluctuating intensity of fluorescence emitted by the molecules or other particles 
in at least one measurement volume of a non-uniform spatial brightness profile by 
measuring numbers of photon counts in primary time intervals by a single or more photon 
detectors, 

b) determining at least one distribution function of numbers of photon counts, 'J^n), from the 

measured numbers of photon counts, 

c) determining physical quantities characteristic to said particles by fitting the experimentally 
determined distribution function of numbers of photon counts, 



■ 

wherein the fitting procedure involves calculation of a theoretical distribution function of the 
number of photon counts P(ik) through its generating function, defined as 



2. A method according to claim 1 wherein the primary tune intervals are consecutive intervals 
of equal width, 

3. A method according to one of claims 1 and 2 wherein in step b) numbers of photon counts 

{n t } subject to determination of a distribution function S(n) are derived from numbers of 

photon counts in primary time intervals {Nj} by addition of numbers of photon counts from 
primary time intervals according to a predetermined rule. 



4. A method according to claim 3 wherein in step b) numbers of photon counts subject to 
determination of a distribution function ^(n) are calculated from the numbers of photon 

counts in primary time intervals {Nj} according to the rule n, =^N Vfi+ft , where M is an 

integer number expressing how many times the time interval in which {nj} is determined is 
longer than the primary time interval. 



5 A method according to claim 3 wherein in step b) numbers of photon counts <*> axe 
time intervals axe separated by a time delay. 

' T ' ' ' ■ * * 1 

* ' ' ^ ' " ■ . '■*.'' " ' ' 

6 A meAod according to claim 5 herein in step b) numbers of photon counts {*) subject to 
: det erh^ °* *** ^ nurabers of^photon 

counts in primary time intervals m *> > ^ =g(N^ + N^^ 
where M and I are positive integer numbers. : / ^ ^ 

,. A nKthodKcording to one 6f the cWms 3 lo 6 hereto in step b) a set ° f aisoa »" i< " > 
(tactions is determined .ccorimg to a set of different rules. »d set of distribution 

V functions being fitted jointly in step c). - 

• , •■- n w v p> _,; n - lft steo cl a set of distribution functions with 

8. A method according to ciaam 7 wherein in step cj a . 

different values of M and/or I are fitted jointly. 

9. A method according to one of the claims 1 to 8 wherein at least one of the physical 
quantities of step c) is concentration of particles. 

10. A method according to one of the claims 1 to 9 wherein at least one of the physical 
quantities of step c) is specific brightness of particles. \ ; - ; V 



IL A method according to o ne of the claims 1 to .10, wherein.at least one of the physical 
quantities of step c) is diffusion coefficient. 

V A method according to one of the claims 1 to 11 wherein the generating function is 
calculated using the expression G<^exp[J«*^^ 

the density of particles with specific brightness *, T is the length of the counting interval, 
and B(r) is the spatial brightness profile as a function of coordinates. 
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13. A method according to one of the claims I to 12 wherein the argument of the generating 
function is selected in the foim § = e~* and a fast Fourier transform algorithm is used in 

calculation of the theoretical distribution of the number of photon counts out of its 
generating function- 

14. A method according to one of the claims 1 to 13 wherein in step c) when calculating the 
theoretical distribution <P(n), the spatial brightness profile is modelled by a mathematical 
relationship between volume and spatial brightness, 




15. A method according to claim 14 wherein in step c) when calculating the theoretical 
distribution P(n\ the spatial brightness profile is modelled by the following expression: 

dV 13 

= a,x + a 2 x + a % x , where dV denotes a volume element, x denotes logarithm of the 

relative spatial brightness, and a u a 2 and « 3 are empirically estimated parameters. 




16. A method according to any of the claims 1 to 15 wherein in step a) a confocal microscope 
is used for monitoring the intensity of fluorescence. 

17. Use of a confocal apparatus for performing the method according to any of the. claims 1 to 
16 comprising: 

a radiation source (1 2) for providing excitation T&diation (14), 

an objective (22) for focussing the excitation radiation (14) into a measurement volume 
(26), 

a detector (42) for detecting emission radiation (30) that stems from the measurement 
volume (26), and 

an opaque means (44) positioned in the pathway (32) of the emission radiation (30) or 
excitation radiation (14) for erasing the central pan of the emission radiation (30) or 
excitation radiation (14), 
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A method for characterizing 
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fluorescent molecules or other particles in samples comprising 



a 



^ ***** <— * -^*h**** 

ta W M on. m^cm vohrne of . non-uniforn, *U ongh^s proffle .* 
curing ambers of photonooun,™ *- * ««•*«» 

detectors, 

b> *«**« » 1«* one dtobmion tat*, of numbe* of photon «.»«. #»), fron,the 

measured numbers of photon counts, 
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